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ABSTRACT 

Taking into account noise from intrinsic ellipticities of source galaxies, in this 
paper, we study the peak statistics in weak-lensing convergence maps around 
clusters of galaxies and beyond. We emphasize how the noise peak statistics is 
affected by the density distribution of nearby clusters, and also how cluster-peak 
signals are changed by the existence of noise. These are the important aspects 
to be understood thoroughly in weak-lensing analyses for individual clusters as 
well as in cosmological applications of weak-lensing cluster statistics. We adopt 
Gaussian smoothing with the smoothing scale 6g = 0.5 arcmin in our analyses. 
It is found that the noise peak distribution near a cluster of galaxies depends 
sensitively on the density profile of the cluster. For a cored isothermal cluster 
with the core radius Rc, the inner region with R < R^. appears noisy containing 
on average ~ 2.4 peaks with i/ > 5 for R^. = 1.7 arcmin and the true peak height 
of the cluster u = 5.6, where u denotes the convergence signal to noise ratio. For 
a NFW cluster of the same mass and the same central u, the average number 
of peaks with u > 5 within R < R^ is ^ 1.6. Thus a high peak corresponding 
to the main cluster can be identified more cleanly in the NFW case. In the 
outer region with Rc < R < 5Rc, the number of high noise peaks is considerably 
enhanced in comparison with that of the pure noise case without the nearby 
cluster. For z/ > 4, depending on the treatment of the mass-sheet degeneracy in 
weak-lensing analyses, the enhancement factor / is in the range of ~ 5 to ~ 55 
for both clusters as their outer density profiles are similar. The properties of the 
main-cluster-peak identified in convergence maps are also affected significantly 
by the presence of noise. Scatters as well as a systematic shift for the peak 
height show up. The height distribution is peaked at u ~ 6.6, rather than at 
u = 5.6, corresponding to a shift of Au ~ 1, for the isothermal cluster. For the 
NFW cluster, Au ~ 0.8. The existence of noise also causes a location offset for 
the weak-lensing identified main-cluster-peak with respect to the true center of 
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the cluster. The offset distribution is very broad and extends to i? ~ i?c for 
the isothermal case. For the NFW cluster, it is relatively narrow and peaked 
at i? ~ 0.2i?c- We also analyze NFW clusters of different concentrations. It is 
found that the more centrally concentrated the mass distribution of a cluster is, 
the less its weak-lensing signal is affected by noise. Incorporating these important 
effects and the mass function of NFW dark matter halos, we further present a 
model calculating the statistical abundances of total convergence peaks, true and 
false ones, over a large field beyond individual clusters. The results are in good 
agreement with those from numerical simulations. The model then allows us 
to probe cosmologies with the convergence peaks directly without the need of 
expensive follow-up observations to differentiate true and false peaks. 

Subject headings: cosmology: theory — dark matter — galaxy: cluster — general 
— gravitational lensing — large-scale structure of universe 

1. Introduction 

Arising from the light deflection by the inhomogeneous mass distribution of the universe, 
gravitational weak-lensing effects have been applied to map out, individually or statistically, 
the mass distribution of galaxies, clusters of galaxies, and large-scale structures (e.g., Bartel- 
mann & Schneider 2001; Mandclbaum ct al. 2008; Li ct al. 2009; Clowe et al. 2006; Corless 
et al. 2009; Fu et al. 2008). Deep observations also reveal the time evolution of the structure 
formation (Massey et al. 2007). Cosmic shear analyses have been used to constrain differ- 
ent cosmological parameters (e.g., Benjamin et al. 2007; Kilbinger 2009; Li et al. 2009). 
Future wide and deep weak-lensing surveys will give rise to observational data with greatly 
improved quantity and quality. Therefore the derived cosmological information is expected 
to be much more accurate than that of today, which will make the weak-lensing effect one 
of the most important probes in cosmology, notably, in understanding the nature of dark 
energy (e.g., Albrecht et al. 2006). However, it has been realized that weak-lensing effects 
are not as clean as they appear to be. Being very weak, they can be affected and contami- 
nated easily by different observational and physical effects. Considerable efforts have been 
paid to understand these effects and their influences on cosmological studies (e.g.. Tang & 
Fan 2005; Heymans et al. 2006; Ma et al. 2006; Huterer et al. 2006; Bridle & King 2007; 
Fan 2007; Amara & Refregier 2008; Sun et al. 2009). 

Among others, intrinsic ellipticities of source galaxies are known to dominate over their 
weak-lensing induced shear signals by a large factor. The environmental dependence of 
galaxy formation can lead to intrinsic-intrinsic and shear-intrinsic correlations, which can 
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contaminate directly the two-point weak-lensing shear correlations at percent to tens of per- 
cent level (e.g., Hirata & Seljak 2004; Mandelbaum et al. 2006; Fan 2007; Zhang 2008; 
Joachimi & Schneider 2009). Moreover, in reconstructing the foreground mass distribution 
from shear measurements of background galaxies, their intrinsic ellipticities leave consider- 
able noise even if they are completely random (e.g., van Waerbeke 2000; Hamana et al. 2004; 
White et al. 2002). Understanding fully the statistical properties of the noise field and the 
mutual influence between noise and the true lensing effects are therefore crucially impor- 
tant in studies aiming to reveal the mass distribution of large-scale structures. Particularly, 
chance alignments of intrinsic ellipticities can form false peaks in a lensing convergence map, 
and mimic 'dark clumps', namely structures with large mass to light ratios (e.g., Erben et 
al. 2000; Linden et al. 2006; Gavazzi & Soucail 2007; Schirmer et al. 2007; Fan 2007). Fur- 
thermore, noise can affect true lensing signals, causing large scatters and possibly systematic 
biases to them. These can lead to considerable uncertainties in analyzing the mass distribu- 
tion of clusters of galaxies, in investigating the possible existence of 'dark clumps', and in 
deriving constraints on cosmological parameters from weak-lensing cluster abundances. 

Different statistical methods have been developed to suppress such noise in weak-lensing 
mass reconstructions, and the smoothing is one of the simplest ways for doing so. Given the 
rms of intrinsic ellipticities, the amplitude of the smoothed noise depends on the smoothing 
scale and on the surface number density of source galaxies thus the observational depth 
(e.g., van Waerbeke 2000). Because of the central hmit theorem, it is expected that the 
statistics of the noise field after smoothing should be close to Gaussian if a large enough 
number of galaxies are included within the smoothing window. By employing a Gaussian 
smoothing function, van Waerbeke (2000) shows that the noise field of the convergence is 
indeed approximately Gaussian. He further analyzes its properties, such as the number and 
the shape distribution of the noise peaks. Fan (2007) extends the study by including the 
intrinsic alignment of background galaxies in the analyses, and shows that the number of 
false peaks are enhanced depending on the strength of the intrinsic alignment. Schneider 
(1996) consider the effects of noise on the height of the peaks associated with true mass 
concentrations in aperture- mass maps by adding a Gaussian scatter to the peak height. 
From numerical simulations, Hamana et al. (2004) analyze noise effects on true peaks, and 
present a phenomenological correction to the abundance of weak-lensing detected peaks. In 
this paper, we put forward a theoretical model to study the mutual influence of noise and 
true peaks. Specifically, we consider the lensing effects of a cluster of galaxies taking into 
account noise from intrinsic ellipticities of source galaxies. We investigate the effect of the 
density profile of the cluster on the distribution of noise peaks, and conversely, how noise can 
change the true peak signal of the cluster. Based on these analyses and the mass function 
of dark matter halos, we develop a model to compute statistically weak-lensing convergence 
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peak abundances over a large field, including both the peaks associated with clusters of 
galaxies and the false ones from noise. Therefore we can use convergence peaks directly 
as cosmological probes without the need to sort out true and false peaks. Moreover, the 
distribution of noise peaks depends on the density profile of true clusters, and thus contains 
additional cosmological information, which is naturally included in our model. 

The rest of the paper is organized as follows. In §2, we present our theoretical framework. 
The results are shown in §3 with §3.1 for individual cluster cases and §3.2 for statistical results 
on a large field. §4 contains summaries and discussions. 



2. PeELk statistics in noisy convergence maps around clusters of galsLxies 

With the intrinsic ones taken into account, the observed ellipticities of source galaxies 
can be written in the weak lensing regime as (e.g., Bartelmann & Schneider 2001) 

where 7 is the complex lensing shear due to the foreground mass distributions, and e*^*-*) 
and e^^^ represent the observed and the intrinsic complex ellipticities of source galaxies, 
respectively. The corresponding smoothed quantities are 

1 

where E*^*^^ and T are the smoothed e*^°^ and 7, respectively, W(6) is the smoothing function, 
and Ug and Ng are, respectively, the surface number density and the total number of source 
galaxies in the field. 

The lensing convergence k can be obtained from the shear through their following rela- 
tion in the Fourier space 

Hk) ^ Ca{k)%{k), (3) 

— * 

where the summation over a — (1,2) is implied, and Cq = [cos(2(/7), sin(2(/7)] with k — 
k {cos (p, sirup) (Kaiser & Squires 1993). Then the real-space convergence reconstructed from 
the observed ellipticities is 

Knie) ^ Jdk e-*^"'c„(fc)e(p)(fc) = k{9) + n{e), (4) 

where k and n are the lensing convergence and the noise from the intrinsic ellipticities, 
respectively. After smoothing, we have 

KN{0)^K{e) + N{0) (5) 
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with K and N the corresponding smoothed k and n. It is therefore clear that in order to 
extract the true K from the noisy Kjv, we must understand the influence of N . 

In this paper, we mainly consider foreground clusters of galaxies as lens objects. Con- 

— * 

sidering an individual cluster, we take K{9) to be the smoothed convergence of the cluster, 
which is determined by its assumed known density profile. The noise is regarded as a 

— * 

random field. The total K]\f{6) is therefore also a random field with its statistics deter- 
mined by the noise field N{6) and modulated by the cluster K{9). Concerning peaks in 
the noisy convergence Kj<q{d) around an individual cluster, there should be a dominant one 
corresponding to the main cluster. There are also false peaks arising from A^. It is expected 
while the properties of the dominant peak in Kj\f[9) should be closely related to those of 
true K{9), they are inevitably affected by the existence of N. On the other hand, the noise 

— * — * 

peaks in Kn{9) should also be modulated by the cluster K{9) in comparison with those 
in pure noise field N. It should be noted that the peak statistics depend not only on the 
value of Kn, but also on its profile through the first and second derivatives of K^. Thus 
K and need to 'collaborate' to give rise to peaks in Kn. Therefore the peak properties 
in Ki\i differ from those in iV or in in ways that are more complex than just changing 
the peak height by a value of K or A^. Our main focus in the paper is to understand such 
mutual infiuences between K and A^ in terms of the peak properties in Kj^, and further their 
effects on relevant cosmological studies. The formulations presented in the following are for 

— * 

noisy Kn(9) around a cluster with a known density distribution. These set the theoretical 
framework for analyses presented in §3.2 on large-scale statistics beyond individual clusters 
incorporating the mass function of dark matter halos. 

It has been shown that the smoothed noise field N is approximately Gaussian because 
of the central hmit theorem (van Waerbeke 2000; Fan 2007). For an individual cluster, its 
convergence K is regarded as a known quantity. Thus the statistics of the field is also 
Gaussian and follows that of N with the modifications due to the presence of K from the 
cluster. For the random Gaussian field Kn, its peak statistics involves random variables 
Kn, diKjN (denoted as Klf) and dijKjy (denoted as K]^) = 1, 2). Their joint probability 
distribution function can be readily written out as follows (e.g., Bardeen et al. 1986; Bond 
& Efstathiou 1987) 

-2 



p{Kn,K]^,K^^,K'^,KIj,KI) dKN dK]^ dK^^ dK^^ dK], dKj, 

1 r {{Kn - K)/ao + MiKjl " K'') + {K]^ - K^')]/a,} 



2 . 



exp i rrr I X 



[27r(l-7^)ao]V2 2(1-7^) 
1 r [{K},' - K'') - {K§ - K'')]' {K]^-K''f {K^^ - K^'f 



27r(7| { 2(tI (tI (jI i ^ 

12 7>'12\2^ 
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— ^ exp 



(Kk - K^f {Kl - K 



2\2 



dKn dKll dK]^ dK^^ dK^ dK^, 



(6) 



where all the quantities related to Kn are attached a subscript N, while those without the 
subscript refer to the ones from K. The ctq, cti and (72 are the moments of the noise field N, 
which are given by, with N{k) being the Fourier transform of N (e.g., van Waerbeke 2000), 



dk k'' < \N{k)\' > . 



(7) 



The quantity 7jv is defined as '-/n = o"^/(o"oO"2). By diagonalizing the second derivative tensor 



—K]^ to obtain the eigen values of Aati and AAr2 and by defining 



02 



(Aati — AiV2) 



(8) 



we have 



K 



11 

N 



a2XN 



(l+2eivcos2^Ar), K 



22 
N 



(72 a; AT 



1— 2ejv cos26'jv), K = —a2Xjqejqsm26i^i, 

(9) 

where 6^ is the rotation angle in the range of [0,7r]. We further define vn = K^/ao, then 
Eq. (6) can be simplified to 

p{uN,XN,e]\f,9]\f, KIj, Kjj) dvN dx^ de^ dd^ dK\i dK^f — 
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{K\, - Kr {Kl - K 



2\2 



dKidKiy 



(10) 



Considering peak statistics, the average number density of maxima given pn = i^o,Xn ~ 
xq, cn = Co, and 9n = 9q is (Bond & Efstathiou 1987) 

npeak{j^o, xo, eo, 6*0) = < S{iyN - J^o)^{xn - Xo)S{eN - eo)S{9N - 9o) 

5(Kl,)5iKl,){al/A)x%(l - 4e^)e(l - 26^^)0(6^^) >, (11) 



where the average is calculated with the probability function given by Eq. (10), and the 
step functions 6(1 — 2eAr) and 0(eAr) are required to satisfy the conditions of maxima with 
Aati > Ajv2 > 0. Thus the number density of peaks given i/jv = i^o and — Xq is 



npeakV^O, ^Oj — ©Xp — ^ 2 ^ ^ CXp 

L Z 0*9 (7*9 (T9 -I L 
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cxp ( — Sxoejv sin 2^jv 
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0-2 



where 9^ — 2(7f/(T|. We reorganize the above equation and obtain the following 



ripeaki'^O: Xo) ^ CXp 



+ iK'yi( 1 1-1 r 1/ ir\2i 
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(13) 
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Further we have 
ripeakii^o) = exp 



(Xl)2+(X2)2 



exp 
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It is seen clearly that the number of peaks near a cluster of galaxies depends on the 
density profile of the cluster through K\ and K'^^ = 1,2). Assuming a spherical 
mass distribution for a cluster of galaxies, we can calculate Upeaki^o) at different radii along 
a specific direction where K^'^ = in our chosen coordinate system. Then we have 

(^11 _ ^22)2. 



F(xn) 



exp 

1/2 



X 



dcN 8(x^eAr)xjv(l - 4ejv) exp(-4x^e^) x 



d9 



N 



TT 



exp 



— Axn^n cos2^iv 



(^11 _ ^22) 



0-2 



(16) 



Without the foreground cluster, the function F{xn) can be written out analytically as 



F{xn) ^ x%-1 + exp{-x%) (Bond & Efstathiou 1987). With a non-zero K^^ 
find a fitting formula for F{xn) with 



K 



22 



we 



F{xn) 



X 



N 



1 + exp( 



-X 



N. 



+ CUeM-^%/9{Ci)]-l}, 



(17) 



where Ci = (K — K )/a2 and g{Ci) is a fitted value depending on Ci. In Figure 1, we 
show two sets of results for F{xiy) — [x% — 1 + exp(— a;^)] with Ci = 0.2 (upper set) and 
Ci = 0.5 (lower set), respectively. For each set, both the numerical result (solid line) and the 
one calculated from the fitting formula (dashed fine) are shown. For Ci = 0.2 and Ci = 0.5, 
we have g{Ci) — 3.3 and g{Ci) — 3.5, respectively. It is seen that Eq.(17) provides an 
excellent approximation to F{xn), which will be adopted for further calculations. 

Given a spherical cluster profile, we can obtain its convergence k, and further the 
smoothed K with a smoothing function W through 



K{x)^ 

The relevant derivatives are 

K' = diK{x) 

and 



dx' W{x' — x)k{x'] 



dx diW{x' — x)k{x'), 



K'^ = dijK{x) = J dx' dijW{x' - x)k{x'). 

In our analyese, we use the Gaussian smoothing function with 

1 



W{x) 



\x\'^ 



exp 



Oh 



(18) 

(19) 
(20) 

(21) 
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where Oq is the smoothing scale. We calculate K, and K'^^ along the xi (i.e, X2 — 0) 
direction. Then the relevant quantities become 



01 



1-1 -^Di 



(22) 
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dec dx 2 
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X exp 
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(25) 



and ii-^ = and K'' = 



We consider two types of most commonly used cluster profiles, namely, the isothermal 
profile and the NFW profile (Navarro, Prenk & White 1996). For the isothermal one, we 
adopt the non-singular isothermal sphere model with its lensing potential having the form 
(van Waerbeke 2000) 



0(Xi,X2) = A^A 1 + 



[Xi - XioY + {^2 - X20) 



(26) 



where Rc is the core radius, {xiq,X2o) is the central position of the cluster, and 

= {47ra^Jc^){DoiDijD,,)Rl. (27) 

Here Doi, Dis and Dos are the angular diameter distances between the observer and the lens, 
between the lens and the source, and between the observer and the source, respectively, and 
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(7c represents the velocity dispersion of the isothermal cluster. We choose {xio, X20) — (0, 0), 
then the convergence k of the cluster is 



Ko 2 + {xl + xl)/Rl 

K,{Xi,X2) - 



2 [l + (x2 + a:2)/i?2]3/2' 

where kq is related to and Rc. 

For the NFW profile, we have (Navarro, Frenk & White 1996) 

Ps 



(28) 



p(r) = 



r/rs(l + r/r5)2' 



(29) 



where ps and are the characteristic mass density and scale of a dark matter halo. The 
corresponding convergence is k,{R) — {2psrs/T,c)f{R/rs) with 



1 - {R^f 



+ 



1 - (R^) 



3/2 



1 

3 



1 - {Rsf 



(i?^) -1 



3/2 



arccosh(l/^*) for < 1 

for = 1 
arccos(l/^') for > 1 (30) 



where i?* = \/ x1 + x1/rs. Given the redshift of the lens, i?** can also be expressed in angular 
coordinates with W = ^/(xi/A)^ + (2:2/^)^/ ijs/Di), so for the isothermal case of Eq. (28). 
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In Figure 2, we show the smoothed quantities Kja^^ jo\^ {K^^ + K'^'^)/ (J2. and {K 
K'^'^)/a2 for the two profiles, respectively. The specific parameters are Oq = 0.5 arcmin, 
Rc = 1.7 arcmin and kq/cxq = 6 for the isothermal model (solid lines). For the NFW model 
(dashed lines), we adjust Rg and the central amplitude of k,{R) so that it has the same mass 
and the same smoothed central K as those of the isothermal cluster. In this case, Rg — 5Rc. 
Note that the horizontal axes in the plots are R — R/Rc- It is seen that the two profiles 
are different, especially in the inner region. Thus we expect to see differences of the peak 
statistics in the two cases. 



3. Results 

This section contains two parts. In §3.1, we show the results around individual clusters. 
The isothermal and the NFW clusters are analyzed. The field extends to 5Rc for all the 
considered cases. The analyses shown here are relevant to weak lensing studies targeted at 
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individual clusters. In §3.2, the statistical results over a large field are shown. Specifically, 
we present a theoretical model based on the mass function and the NFW density profile 
of dark matter halos to predict the numbers of high peaks in the weak lensing convergence 
map over a large field. Our model takes into account the effects of noise on the peak 
heights of true clusters and the enhancement of noise peaks near clusters of galaxies. Given 
a cosmology, our model can then predict the number distribution of peaks with different 
heights. These cosmology-dependent quantities can be used as very promising and efficient 
probes to constrain cosmologies. The advantages of these probes are two folds. Firstly, 
because we concern the total number of peaks, we do not need to differentiate true peaks 
and false peaks, which often requires extensive follow-up observations. Secondly, the redshift 
information of the peaks is not necessary since we use the peak height distribution rather 
than the redshift distribution of peaks as cosmological probes. 



3.1. Individual cluster fields 

Here we discuss the peak statistics around individual clusters of galaxies. Two aspects 
are addressed. One is how the distribution of noise peaks is influenced by the density profile 
of a cluster. The other is how the properties of the peak corresponding to the true cluster, 
such as its height and position, are affected by the presence of noise. In our calculations, 
we assume that a considered cluster is located at the center, and the field extends to 5i?c- 
Our main results are presented for the two clusters shown in Figure 2. In addition, we also 
analyze different NFW clusters of the same mass but with different concentrations to further 
demonstrate the dependence of the peak statistics on cluster profiles. Besides theoretical 
analyses, we perform Monte Carlo simulations with different realizations of the noise field 
added to the considered cluster re-map. The theoretical results are then compared with the 
results averaged over 100 simulations for each of the clusters. 

For a visual impression, we first present in Figure 3 a typical set of maps from our Monte 
Carlo simulations for the two clusters shown in Figure 2. Figure 3a and Figure 3b are for the 
isothermal cluster, and the NFW cluster, respectively. It is seen that in the regions outside 
Rc, there is almost a one-to-one correspondence between the peaks in the two maps. This 

is because in these outer regions, the density profiles of the two clusters are relatively flat 
with small values of K^/ai, K^^ + A'^^/o"2 and K^^ —K'^'^joi. Then the peaks are essentially 
the ones existed in the pure noise field but with their heights modified by the term K/cjq. 
However, it should be noted that K/ao depends on the radius i?, and thus the peak statistics 
over a range of R cannot be described by an overall constant shift in peak height from that 
of N. Within R^ the peaks look very different in the two maps, which clearly shows the 
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effects of the density profile (K^/ai, K^^ + K'^'^/a2 and K^^ — K'^'^ / in addition to tfie term 
K/aQ. For tlie NFW cluster, there is a single high peak which corresponds clearly to the 
peak signal of the main cluster. For the isothermal cluster, however, two prominent peaks 
with comparable heights occur. The positions of the peaks are also shifted from the center. 
Such peak configurations due to noise can cause complications in weak lensing studies of the 
structure of clusters. 

To be quantitative, we present in Figure 4 the radial distributions of the number of 
peaks above a certain threshold around the two clusters. In each panel, the lines are the 
results from our theoretical calculations, and the symbols arc the average results from the 
100 Monte Carlo simulations for each case. The dash-dotted line and the '+' symbols are 
for the NFW cluster. The sohd fine and the diamonds are for the isothermal cluster. The 
vertical axes are for the number of peaks within rings of width = 0.1, where R — R/Rc- 
Different scales of the horizontal and vertical axes in different panels should be noted. It is 
seen that the results from the theoretical model and from the simulations agree very well. 
The peak distributions for the two cluster profiles are distinctively different, especially in the 
inner region and for high peaks with z/ > 5. For the NFW cluster, high peaks are mainly 
found near the central region of the cluster, and the total number of high peaks is smaller 
than that of isothermal cluster. In the region with R < Rc, the average numbers of peaks 
for i/ > 3,4,5 and 6 are approximately 3, 2.6, 1.7 and 0.8, respectively, for the NFW case. 
For the isothermal case, the corresponding numbers are 3.5, 3.2, 2.4 and 1.1. Therefore, for 
the NFW cluster, the inner region is relatively clean and often contains a single and well 
defined high peak as seen in Figure 3b. On the other hand, for the isothermal cluster, there 
are on average 2 peaks with z/ > 5 in the inner region, of which, one corresponds to the real 
cluster and the other is generated due to noise (Figure 3a). 

We further consider the inner {R < Rc) and outer {Rc < R < 5Rc) regions separately. In 
the outer regions, we analyze the occurrence probability of high noise peaks, which is relevant 
to the dark clump problem. We emphasize, with respect to the case without clusters, the 
enhancement due to the cluster density distribution. Figure 5a shows the distribution of the 
cumulative number of peaks above the threshold Uf in the region 1 < i? < 5 (7? = R/Rc) 
and Figure 5b presents the corresponding enhancement factor / = Npcak{^t) / ^rani^t) i where 
Npcaki^t) and N^ani^t) denote, respectively, the cumulative numbers of peaks in the cases 
with and without a cluster in the center. In Figure 5a, the solid, dash-dotted and dashed lines 
are for the isothermal cluster, the NFW cluster, and the case without a cluster, respectively. 
The diamond symbols and the plus symbols are for the results averaged over 100 simulations 
of the isothermal and NFW clusters, respectively. In Figure 5b, the solid and dash-dotted 
lines are for the corresponding enhancement factors of the two clusters. It is seen that the 
results are similar for the two clusters because their density profiles are similar in the outer 
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regions as shown in Figure 2. However, differences can still be seen for high peaks. For 

Vt > 5, the enhancement factor is larger for the isothermal cluster than that for the NFW 
cluster. This is due to the fact that in the outer regions, these high peaks are mainly found 
in 1 < i? < 2, where K/(Tq is higher for the isothermal cluster (Figure 2a). For Ut = 3, 4, 5, 
6, / ~ 9.8, 55, 595 and 10000 for the isothermal cluster. The corresponding / for the NFW 
cluster are 10.7, 55, 437, and 4671. Therefore depending on the density profile of a cluster, 
the occurrence probability of high noise peaks around it can be greatly boosted. The effects 
are larger for higher peaks. 

As discussed earlier, in the outer region, the effect of the density profile of a cluster 
on noise peak statistics comes in mainly through the term K/ao- It increases the heights 
of the peaks originally with — K/ao < v < Vt to he above the threshold fj. Because 
of the approximately exponential dependence of Npeaki^^) on u, the enhancement ai Ut > 3 
changes sensitively with the values of K/ao and Ut. This enhancement in outer regions can 
be analogous to the peak-background split scenario in explaining the biased dark matter 
halo distributions with respect to the dark matter mass density distribution (Bardeen et al. 
1986; Mo & White 1996). It is noted, however, that K/ao changes with radius. Thus the 
local enhancement depends on the radius, and the overall enhancement factor / over the 
outer region with 1 < < 5 cannot be modeled by a constant shift of with respect to 
that of A^. When the observational field is comparable to the size of a cluster, we should 
also be aware that the mass-sheet degeneracy in lensing analyses brings some arbitrariness 
in the reconstructed mass distribution, which in turn affects the noise peak statistics within 
the field. In the weak lensing regime, this arbitrariness is mainly reflected by an additive 
constant to k. In real studies, different people may choose different constants. It can be 
chosen so that the averaged < k, >= in the observed field (e.g., Erben et al. 2000), or 
to make the k averaged over the outer-most ring of the observed field be zero (e.g.. Kaiser 
& Squires 2003). It is also attempted to recover the true mass distribution of a cluster by 
adjusting the constant to make the reconstructed mass distribution be consistent with the 
distribution determined by some other observations (e.g., Umetsu & Futamase 2000). In 
Figure 6a and Figure 6b, we show / for the isothermal and NFW clusters, respectively, by 
applying an additive constant to the k field to have < k >— over the area of < 5 
(dashed), to have < k >= at i? = 5 (dash-dotted). The solid lines are for the true mass 
distribution. It is seen that different choice of the constant results significantly different 
enhancements. However, even in the case of < k >= over the field within R = 5, which 
gives the lowest enhancement factor, / ~ 10 and 6.2 at i^t = 4.5 for the isothermal and NFW 
clusters, respectively. 

The strong association of noise peaks with their nearby clusters can have significant 
effects on cosmological studies with weak lensing effects. It has been argued that the 'dark 
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clump' found near cluster A1942 (e.g., Erben et al. 2000) is unlikely to be a noise peak 
resulting from the chance alignment of source galaxies. This is based on the peak statistics 
without considering the cluster environment, which gives the average number of noise peaks 
with V > 4.5 in the considered area of < R < 5Rc to be about 0.01 for 9g = 0.5 arcmin. 
From our results shown here, the number of noise peaks can be increased by a factor of 
/ ~ 10 if the density profile of the cluster is approximately isothermal. This leads to an 
order of magnitude increase in the probabihty that the seen 'dark clump' near A1942 is a false 
one without corresponding to real mass concentrations. Thus in order to properly quantify 
the statistical significance of a 'dark clump' near a cluster, the cluster mass distribution has 
to be taken into account. Specifically, in observational weak lensing studies, such as the one 
for A1942, one can first make an approximate modeling of the mass distribution of the main 
cluster from the weak lensing measurements. Then our model can be applied to study the 
occurrence probability of noise peaks in its surrounding area. Such a probability is more 
appropriate than the one from the blind estimate without considering the existence of the 
cluster itself. It should also be emphasized that this probability depends sensitively on the 
density profile of the cluster. Moreover, because of its strong dependence on the density 
profile of the cluster, the spatial distribution of noise peaks can mimic, to a certain extent, 
the true substructure distribution within the cluster. Thus the noise peaks can contaminate 
the cosmological information inferred from weak lensing substructure analyses. 

We now turn to the inner region with R < Rc or R < 1. The highest peak within this 
region is assumed to correspond to the signal from the cluster itself. We then analyze how 
the height and position of this main-cluster peak are affected by the existence of noise. 

In the case without noise, the expected weak lensing peak heights should be the same 
for a sample of clusters with the same mass distribution and at the same redshift. With 
noise included, however, their measured peak heights do not have the same value anymore. 
Instead, they follow a probability distribution. Based on our analyses shown in the previous 
section, the probability can be written as 



where Upeaki^) is given by Eq.(15) and the central values of K, K^, K'^, K^^ and K^"^ at i? = 
are used in calculating the probability. In Figure 7, we show p{v) for the isothermal cluster 
(upper panel) and the NFW cluster (lower panel), respectively. In each panel, the solid line 
is for the theoretical result from Eq. (31), and the histograms are for the result from the 100 
Monte Carlo simulations. For each simulation, we identify the highest peak within the region 
^ < 1 to be the one corresponding to the main cluster. The true cluster mass distribution 
is considered without additional treatments for the mass-sheet degeneracy. Firstly, it is seen 
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that the sohd hne is in very good agreement with the results from simulations, showing that 
Eq.(31) can properly model the effects of noise on the measured height of the main-cluster 
peak. Secondly, it is noted that the average value of the measured peak height is shifted 
toward larger v with respect to the true height of the cluster (indicated by the vertical 
dotted line in each panel). The specific value of the shift depends sensitively on the density 
profile of the cluster. For the isothermal cluster, the shift is about Az/ ~ 1.1 from v = 5.6 
of the true peak height to u — 6.7, whereas for the NFW cluster, the shift is smaller with 
Au ~ 0.8. This systematic shift due to the existence of noise is closely related to the fact 
that the differential peak number distribution of a pure Gaussian noise field (considering 
only maxima peaks) is peaked at ~ 1, rather than at u — 0. Such a positive shift can 
increase the number of weak-lensing cluster detections, and thus affect the corresponding 
cosmological applications considerably. In the next subsection, we present a detailed model 
to calculate quantitatively the expected number of weak-lensing convergence peaks, true 
or false, from a large-scale weak lensing survey, taking into account such a shift and the 
enhancement of the noise peaks around clusters of galaxies. This model allows us to be 
able to use directly the high peaks detected in a large-scale weak-lensing convergence map 
as cosmological probes without the need of expensive follow-up observations to differentiate 
true and false peaks. 

The existence of noise also results a spatial offset for the identified main-cluster-peak 
from the true central position of the cluster. This can be understood as follows. Without 
noise, the main-cluster-peak in the weak lensing convergence map locates exactly at the 
cluster center {R — 0) in our considered cases of spherical clusters. At this spatial location, 
the conditions diK = (i = 1, 2) are satisfied. With noise, however, the location of a peak 
should satisfy necessarily the conditions diK^ = di{K + N) = {i = 1,2), rather than 
diK = {i = 1,2). Given a realization of a noise field N and considering a small offset 
ARi {i = 1, 2), the original central peak now shifts to a location which can be estimated by 
{dijK)ARj = —diN = 1,2), where we have used diK^ ~ dijKARj. Thus the offset 
distribution depends on the distribution of diN, which, in the considered Gaussian case, is 
characterized by ai defined in Eq. (7). It also depends sensitively on the density profile 
of the main cluster through dijK. Specifically, the offset distribution can be approximately 
modeled by 

where the peak number density Up^akiAR.v) can be obtained from Eq. (15) by Taylor 
expansion of the quantities of K, and K'^^ with respect to AR. In our considered case. 
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diK{R = 0) = and di{K^^){R = 0) = 0, then to the hnear order of Ai?, we have 

. d^R , A^exp{-[;^n(0)Ai^/.,]n 

In terms of A^, Eq.(33) can be written as 

/Af,N,Af, ,A f. A^exp[-(A^/a)2l 

p(AR)dAR dAR- ^ / \{^^ 34 

^ {0.5a2[l -exp(-a-2)]} ^ ^ 

where = [K^^[0)/a2]{(T2/cTi){Rc)- In our considerations with the Gaussian smoothing 
scale 9g = 0.5 arcmin and Rc = 1.7 arcmin, we have [K^^{0)/(r2] ~ 0.25 and 0.68, and the 
corresponding a ^ 0.58 and 0.22 for the isothermal and NFW cases, respectively. We then 
expect a wider distribution of the offset for the isothermal case than that of the NFW case. 
Figure 8 shows the results of p{AR). The symbols are the simulation results with diamond 
and plus for the isothermal and NFW cases, respectively. The solid and the dash-dotted 
lines are the corresponding results calculated from Eq. (34). Because we use the linear 
expansion in A^ in Eq. (34), we expect some inaccuracies of Eq. (34) when A^ is relatively 
large. On the other hand, we see that Eq. (34) still gives a reasonable description about the 
offset distribution, especially in terms of the location of the maximum probability. For the 
isothermal case, the offset distribution is peaked around AR ^ 0.4, and for the NFW case, 
it is peaked at AR ^ 0.2. The relatively flat density distribution of the cored isothermal 
cluster can therefore lead to an offset of order of Rc for the main-cluster-peak identified from 
weak-lensing analyses. This large offset may have significant effects on the cluster mass and 
the mass profile estimations from weak-lensing observations. For the NFW cluster, the offset 
is smaller. 

The above analyses are mainly done for the two clusters shown in Figure 2, in which, 
the NFW cluster has the same central K and the same mass within R < 5Rc as those of 
the isothermal cluster. This leads to i?^ = 5ii!c, or a concentration parameter c = 1 for the 
NFW cluster. To further demonstrate the dependence of the peak statistics on the density 
profile of a cluster, we also carry out detailed analyses for two more NFW clusters with 
different concentrations and peak heights. The parameters are Rg = l.SRc, and the central 
K = 1.8K{Rs = 5Rc) for one cluster, and Rs = Rc, and the central K = 3K{Rs = 5Rc) for 
the other. These parameters are chosen so that the clusters have about the same mass as 
the NFW cluster with Rs = 5Rc considered above. In Figure 9, we show the results for the 
three NFW clusters. Figure 9a presents the profiles of K/gq. Figure 9b and Figure 9c are 
for the height distributions and spatial offset distributions of the main-cluster peak. In each 
panel, the solid, dashed, and dash-dotted fines are for the cases with Rg — Rc, Rs — 1-SRc 
and Rs — 5Rc, respectively. The histograms in Figure 9b and symbols in Figure 9c are for 
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the results from Monte Carlo simulations. The vertical lines in Figure 9b indicate the true 
cluster peak height (dotted) and the shifted peak height corresponding to the maximum of 
the distribution (dash-dotted) in each case. The dependence of the results on the cluster 
profile is clearly seen. With the same mass, the more centrally concentrated the cluster 
density distribution is, the less the main-cluster-pcak properties are affected by noise. The 
peak height shift due to noise (Figure 9b) is iS.v pa 0.2, 0.4, and 0.8 for the Rs = Rc, 
Rg — l.SRc and Rg = 5Rc, respectively. The corresponding spatial offset distribution is 
peaked at AR ?a 0.03i?c, 0.06Rc and 0.2Rc for the three cases. It is known from Eq. (31) 
and (34) that Au and AR are mainly determined by the central value of \diiK + d22K\ of a 
cluster. A more centrally concentrated density profile leads to a larger value oi\diiK+d22K\, 
and thus smaller Au and AR. 

It is possible to incorporate our model on the height and the position of the main- 
cluster peak [Eq. (31) and Eq. (34)] into observational weak-lensing analyses of cluster 
density profiles so that the noise effects can be properly taken into account. The detailed 
implementations will be explored in our future studies. 

3.2. Peak number statistics in large-scale weak-lensing convergence maps 

Weak-lensing cluster search and subsequently statistical studies are among the impor- 
tant cosmological applications of weak-lensing observations. It is seen from our discussions 
in §3.1 that noise due to intrinsic ellipticities of background galaxies results false peaks in 
the weak-lensing convergence map and the noise peak distributions near clusters of galaxies 
are boosted depending on the density profile of the clusters. It is also shown that the peak 
heights corresponding to clusters of galaxies arc affected by the existence of noise that gener- 
ates large scatters as well as systematic shifts. These effects can contaminate the cosmological 
applications of weak-lensing cluster statistics considerably if they are not taken into account 
carefully. In this subsection, we present a model to calculate the number distribution of 
peaks npeafe(^) over a large area, where npeafe(^) includes both true peaks corresponding to 
clusters of galaxies and the false peaks from intrinsic eUipticities of background galaxies. It 
should be noted that the false peak distribution near clusters of galaxies depends on the 
density profile of clusters, and thus also carries important cosmological information. The 
model then allows us to use directly the peaks detected in the large-scale reconstructed con- 
vergence map from wcak-lcnsing observations as cosmological probes without the need to 
differentiate true or false peaks with follow-up observations. 

In our model, the surface number density of convergence peaks is written as 

npeakWv = nl^^k{v)dv + n'^^^^{v)dv. (35) 
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Here rip^g^f^^u) calculates peaks within virial radii R^ir of dark matter halos, which includes the 
peaks corresponding to halos themselves taking into account the modified peak heights due 
to noise, and the noise peaks with enhanced number distributions near clusters. The term 
^peak^^) noise peaks in the field area away from dark matter halos. Specifically, we 

have 

nUaki^) = J dz'^ J dM n{M, z)f{u, M, z), (36) 

where dV{z) is the cosmological volume element at redshift z, dfl is the sohd angle element, 
n{M, z) is the mass function of dark matter halos given the mass M and redshift z, and 

f{u, M,z)= / dR {2nR) n^eakiy. M, z), (37) 

with n^eakiy^ z) given by Eq. (15). The quantities related to the density profile of dark 
matter halos K\ and K"^^ in Eq. (15) are computed assuming a NEW profile for a halo 
with its characteristic density ps and scale (or the concentration Cyir) determined by M 
and z (Navarro, Prenk & White 1996). 

For n^g,ak{^)i written as 

Keak{^) = ^{^ranil^) [^^ ' j / n{M, z) (vri?^,,)] }, (38) 

where Urani^) is the differential surface number density of noise peaks without foreground 
clusters. 

Wc should emphasize that Eq. (35) does not correspond to a simple and direct summa- 
tion of the number of peaks corresponding to true dark matter halos and the peaks from pure 
noise field. Rather, the two terms in the right-hand side of Eq. (35) represent peaks from 
different environments. The considered area (1 deg^ here) is divided into regions occupied 
by dark matter halos and the field region away from halos. The first term counts for total 
number of peaks in the regions occupied by dark matter halos, which includes both the peaks 
of the halos themselves and the noise peaks within the virial radii of the halos. The halos 
are assumed to be randomly distributed spatially with their average abundance given by the 
mass function n(M, z). The peak abundance within the virial radius of a dark matter halo is 
given by / calculated through Eq. (37) and Eq. (15), which takes into account the effect of 
the halo mass distribution on noise peak statistics and the effect of noise on the peak height 
of the halo itself. To a certain extent, the treatment of the first term is similar to the model 
of halo occupation distribution for galaxies, where the peaks corresponding to halos can be 
analogous to central galaxies and the noise peaks can be analogous to satellite galaxies. The 
second term rip^^f^ calculates the number of peaks in the field region away from dark matter 
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halos, which is the product of of surface number density of pure noise peaks and the spatial 
area not occupied by dark matter halos [cf. Eq. (38)]. 

It should also be noted that to apply Eq. (35), we imphcitly assume that the true 
convergence peaks corresponding to real mass concentrations are linked to individual dark 
matter halos. This is a very reasonable assumption concerning peaks with u > 3. On the 

other hand, for lower peaks, they can arise from the projection effects of large-scale structures 
without being associated with certain primary halos. To analyze these low peaks, it may be 
more appropriate to model the large-scale density field as a Gaussian random field K^ss- 
Then the noise effect on the statistics of these low peaks can be investigated by analyzing 
peaks in the total field Klss + N, which is also a Gaussian random field. In fact, it is possible 
to separate the mass distribution in the universe into highly nonlinear part K which can be 
modeled by dark matter halos and the part associated with large-scale structures that can 
approximately be treated as a Gaussian random field Klss- By studying the peaks in the 
total field K + Klss + N , one can possibly obtain the peak statistics over a wide range of 
peak height, from low to high. We will pursue along this line in our future investigations. 
In this paper, we primarily concern high peaks with > 3, and thus Eq. (35) is used in our 
analyses. 

To validate our modeling of Eq. (35), we also analyze the convergence peak statistics 
from the ray tracing simulations by White & Vale (2004). The simulations we use are for 
the fiat ACDM model with VLm = 0.296, h = 0.7, (Tg = 0.93, and n = 1, where Qm, h, ag and 
n are the dimensionless matter density of the universe, the present Hubble constant in units 
of 100 km/s/Mpc, the rms of the extrapolated linear density fiuctuation over 8 Mpc/i~^ and 
the power index of the initial power spectrum of density fluctuations, respectively. In total, 
16 simulated convergence maps, each with 1024 x 1024 pixels corresponding to a held of 
view of 3 X 3 deg^, are used in our analyses. The source redshift is 2;^ = 1. The details of 
the simulations can be found in White & Vale (2004). For each simulated convergence map, 
we add in noise from intrinsic ellipticities following Hamana et al. (2004). Specifically, a 
random Gaussian noise is added to each pixel with the variance given by 

''^ 2 n,el, 

where is the rms of the intrinsic ellipticity of source galaxies, Ug is the surface number 
density of source galaxies, and Opi^ is the pixel size of the simulated map. We take — 0.4, 
Ug — 30 arcmin"^, and 9pix — 180/1024 arcmin. We then apply a Gaussian smoothing 
to each noisy convergence map. Here we use the smoothing scale 9g — i arcmin, which 
is presumably the most optimal smoothing scale in weak-lensing cluster detections (e.g., 
Hamana et al. 2004). 



(39) 
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In our model calculations, we use the same set of parameters as described above. The 
Sheth-Tormen mass function (1999) and the NFW density profile for dark matter halos are 
applied. 

In Figure 10, we show the surface number density npeafe(i^)/deg^ from our theoretical 
calculations (solid line) in comparison with those from numerical simulations (upper his- 
togram). The simulation results are the results averaged over 16 of 3 x 3 dcg^ convergence 
maps, i.e., we count the peaks over the total area of 16 x 3 x 3 = 144 dcg^, and then divide the 
total peak number by 144 deg^ to obtain the average surface number density of the peaks. 
We see that our model calculations agree with the simulation results very well. In the plot, 
we also show Upeaki^) without considering the noise (dotted line), which is given by 

ripeakii^) = J d^^^ J dM n{M, z) 5[M - M{z, v)\ (40) 

where 5[M — M(2;, i/)] is a 5 function representing the fact that without noise, a dark matter 
halo with given M and z has a determined peak height. The lower histograms are for the 
results from numerical simulations without adding noise in the convergence maps. It is seen 
that this set of histograms is systematically higher than that of the dotted line, indicating 
significant projection effects from large-scale structures. 

Furthermore, to demonstrate the effects of noise on the main-cluster-peak heights, we 
also show (long-dashed line) the results of Upeaki^) ioi peaks associated with true dark matter 
halos but taking into account the change of the main-cluster-peak heights due to noise (cf. 
Figure 7 and Figure 9b). This is given by 

Upeakil^) = J dz^^^ J dM n{M,z) Ppeaki^^, M, z), (41) 

where Ppeak{^, M, z) is given by Eq. (31) with K, K^, and K'^^ determined by M and z. 
Comparing the long-dashed line with the dotted line, we see that at > 3, the number of 
peaks Upeaki^) associated with real dark matter halos increases by a factor of ~ 2. This is 
clearly due to the positive shift of the peak height of a dark matter halo from the effects 
of noise (see Figure 7 and Figure 9b), which boosts certain peaks of relatively small dark 
matter halos initially with u < 3 in the case without noise to u > 3 with the presence 
of noise. This is highly relevant to cosmological studies with the number of true weak- 
lensing peaks, i.e., confirmed peaks by follow-up observations. Without properly taking into 
account the noise effects on the peak heights, Eq. (40) would normally be used to predict the 
abundance of weak-lensing peaks above a certain threshold. By fitting such a prediction with 
the observed abundance of true weak-lensing peaks to constrain cosmological parameters, 
we would obtain a set of biased cosmological parameters. The detailed analyses of such bias 
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effects on cosmological parameters will be presented in our forthcoming paper. It should 
also be pointed out that the differences between the long-dashed and solid lines represent the 
number of noise peaks. These differences are much larger than those calculated from nran(^) 
without considering the enhancement of noise peaks near dark matter halos, indicating the 
importance in taking into account such an enhancement in peak number calculations. 

The results shown in Figure 10 are directly comparable with the ones shown in Figure 
Bl of Hamana et al. (2004). By analyzing the results from numerical simulations, they 
present a fitting correction to the number of peaks taking into account the noise effects 
(the solid line in their Figure Bl). While in qualitative agreements, our results match the 
simulation results better. More importantly, our model [Eq. (35)- (38)] is based on statistical 
analyses of peaks near cluster of galaxies, but not a fitting model to the simulations. The 
excellent agreement between our model prediction (the solid line in Figure 10) and that from 
simulations (the upper histogram) shows that our model can indeed give a reliable prediction 
about the number of convergence peaks in weak-lensing analyses. Thus we can use directly 
the distribution of the total number of peaks, true and false ones, to probe cosmologies. 
This advantage would greatly strengthen the cosmological applications of weak-lensing peak 
statistics. 



4. Discussion 

In this paper, we analyze the peak statistics in weak-lensing convergence maps taking 
into account the mutual influences between the mass distribution of clusters of galaxies and 
noise from intrinsic ellipticities of background galaxies. Concerning the main convergence 
peak corresponding to a cluster, the existence of noise affects both its height and its location, 
resulting scatters as well as a systematically positive shift for its height and an offset distri- 
bution for its position. The effects are sensitive to the density profile of the cluster itself. 
In the two considered cases of the cored isothermal and the NFW clusters with the same 
mass and the same central K, the flat density distribution of the isothermal cluster leads 
to a relatively large shift of the peak height with Au ~ 1 and a wide offset distribution of 
the peak location with AR extending to the core radius Re- For the NFW cluster, its inner 
density profile changes more rapidly and has a relatively large \diiK + ^22-^1 at the center. 
This leads to smaller Az/ and AR, with Au ~ 0.8 and AR peaked at AR ~ 0.2i?c- We fur- 
ther analyze different NFW clusters of the same mass but with different concentrations. It 
is shown that the weak-lensing convergence peak signal of a dark matter halo is affected less 
as its concentration gets higher. We have Au ~ 0.2, 0.4, and 0.8, and AR ~ 0.03i?c, 0.06i?c 
and 0.2i?c, for Rg — Rc, Rs — 1.8i?c and Rg — 5Rc, respectively. 
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For noise peaks, their occurrence probability is boosted near cluster of galaxies. In the 
isothermal case, there are on average about 2.4 peaks with z/ > 5 (note that wc have v = 5.6 
for the considered cluster without noise) in the central region with R < R^. While the 
higher one is usually defined to be associated with the cluster itself, the other is a noise peak 
without any related real structures. This noise peak can be mistakenly identified as a large 
substructure and therefore can complicate the interpretation of the cluster mass distribution 
and further the formation of clusters. For the NFW cluster with Rg — 5Rc, the central 
region is cleaner than that of the isothermal case, and contains on average ~ 1.6 peaks with 
1/ > 5. In outside regions around a cluster with R > Rc, high noise peaks can be confused 
with 'dark clumps', and thus lead to incorrect conclusions regarding the structure formation. 
We find that depending on the density profile of a cluster, the average number of high noise 
peaks in the regions R^ < R < 5Rc can be significantly enhanced. This is locally due to 
the quantity K{R) mainly. But it is important to realize that the i?-dependence, i.e., the 
profile of K{R), plays key roles resulting a net overall enhancement in the considered outside 
region, especially when the mass-sheet degeneracy is taken into consideration. For the cored 
isothermal cluster, the enhancement factor / ranges from ~ 5 to ~ 55 for the number of 
peaks with ly > A, depending on the treatment of the mass-sheet degeneracy. For ly > 5, f 
is in the range of ~ 20 — 600. For the NFW cluster, the enhancement is weaker than that 
of the isothermal cluster for high peaks, with / ~ 3.6 — 55 for > 4, and / ~ 11 — 440 
for 1/ > 5. This enhancement is important in analyzing the statistical significance of 'dark 
clump' occurrence near clusters of galaxies. 

Taking into account the effects of noise on the main-cluster-peak heights and the en- 
hancement of the number of noise peaks near dark matter halos, we propose a model in- 
corporating the mass function of dark matter halos to calculate the statistical abundance of 
convergence peaks over large scales. The model prediction is in good agreement with that 
from numerical simulations. It should be pointed out that because of the mutual effects of 
the mass distribution of dark matter halos and noise, the noise peak abundance also carries 
important cosmological information, especially the information related to the density profile 
of dark matter halos. Our model then provides a theoretical framework to use statistically 
the total number of convergence peaks, including the false ones from noise, as cosmological 
probes. This facilitates the cosmological applications of weak-lensing peak statistics in sav- 
ing the expensive follow-up observations to distinguish true and false peaks. Furthermore, 
the useful cosmological information contained in the noise peak distribution is also included. 

Concerning high peaks with u > 3—4, our model can predict the total convergence peaks 
very well because most of the high peaks are associated with massive virialized dark matter 
halos (and the noise peaks around them). For relatively low peaks, the projection effects 
of large-scale structures can be important, which can change the peak height of clusters 
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as well as produce their own convergence peaks without being associated with virialized 
clusters (e.g., Tang & Fan 2005; Marian et al. 2009, 2010; Maturi et al. 2009). To a certain 
extent, the total weak-lensing convergence can be written as Kjsf = K + Kiss + ^ ■, where 
K and are the respective contributions from virialized dark matter halos and the noise 
as discussed in this paper, and Kiss denotes the large-scale structure contribution beyond 
dark matter halos. On very large scales where density perturbations are approximately linear, 
the effects of large-scale structures can be described by a Gaussian random field reasonably 
well (Maturi et al. 2009). Under this Gaussian assumption, our current model on peak 
distribution near clusters of galaxies can be extended to include the effects of large-scale 
structures by analyzing the sum of the two Gaussian random fields + Klssi which is also 
Gaussian, while treating the cluster i^' as a background. Then Eq. (35) to Eq. (38) can be 
used to study large scale peak statistics, in which the Gaussian field is replaced with the 
total Gaussian field N + Klss- It should be pointed out that physical correlations between K 
and Klss can exist, which may complicate the analyses considerably. The above approach 
essentially divides the structures into two parts, virialized dark matter halos, and linear 
large-scale perturbations described by a Gaussian random field. On scales in between where 
nonlinear effects are significant and the Gaussian approximation of K^ss is not appropriate, 
the modeling on Klss deserves careful investigations. A possible way to study them is to 
extend the mass function of dark matter halos to include such nonlinear but non-virialized 
structures. Detailed investigations on the projection effects of large-scale structures, both 
linear and nonlinear, will be performed as one of our main future tasks. 

Weak lensing effects hold great potentials in cosmological studies. With thorough in- 
vestigations on different systematics, we expect that future weak-lensing observations would 
bring us a greatly improved understanding about the dark side of the universe. 
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10373001, 10533010 and 10773001, and the 973 program No.2007CB815401. Zuhui Fan is 
very grateful for the hospitality of the Institute of Astronomy, University of Cambridge where 
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Fig. 1. — The fitting results for F{x) — [x^ — 1 + exp(— x^)] with Ci = 0.2 (upper set) and 
Ci = 0.5 (lower set), respectively. For each set, the solid line is the numerical result, and 
the dashed line is the fitting result from Eq.(17). The g values are g{Ci = 0.2) = 3.3 and 
g{Ci = 0.5) = 3.5, respectively. 
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Fig. 2.— The profiles of iT/ao, K^/(Ti, {K^^ + K^^)/a2 and (K^^-K^^)/a2 for tlie isothermal 
cluster (solid lines) and the NFW cluster (dash-dotted lines), respectively. The two clusters 
have about the same mass and the same central K/ctq. The Gaussian smoothing scale is 
6g = 0.5 arcmin. The specific parameters for the isothermal cluster is = 1.7 arcmin and 
'to/co = 6 (the smoothed central K/ao ^ 5.6). For the NFW cluster, Rg = 5Rc, and the 
central K/ao ~ 5.6. Here R — R/ Re- The profiles are calculated along xi-axis with X2 — 0. 
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Fig. 3. — The smoothed noisy convergence maps in terms of signal-to-noise ratio K^/aQ for 
the isothermal cluster (Figure 3a, left) and the NFW cluster (Figure 3b, right), respectively. 
The cluster parameters are the same as in Figure 2. The symbols show the positions of the 
peaks, and the circle indicates the radius of Rc- Each map has 1024^ pixels corresponding 
to the total area of 8.5 x 8.5 arcmin^. 
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Fig. 4. — The radial distribution of the cumulative number of peaks with > Ut around 
each cluster. The bin size is 0.1 Rc {R — 0.1). The sohd, dash-dotted lines are the results 
for the isothermal cluster and for the NFW cluster, respectively. The diamond and the plus 
symbols are for the corresponding results averaged over 100 Monte Carlo simulations. The 
cluster parameters are the same as those in Figure 2. 
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Fig. 5. — Fig.Sa (left panel): The cumulative number distribution of noise peaks in the 
area within Rc < R < 5Rc. The solid, dash-dotted and dashed lines are for the isothermal, 
NEW, and the pure noise cases, respectively. The diamond and the plus symbols are for 
the average results from our Monte Carlo simulations of the isothermal cluster and of the 
NFW cluster, respectively. Fig. 5b (right panel): The ratio / = Npeak/Nran in the area within 
Rc < R < 5Rc. The solid and the dash-dotted lines are for the isothermal and NFW clusters, 
respectively. The cluster parameters are the same as those in Figure 2. 




Fig. 6. — The ratio / = Npgak/Nran in the area within R^ < R < 5Rc. Fig. 6a (left panel) is 
for the isothermal cluster, and Fig. 6b (right panel) is for the NFW cluster. In each panel, 
the sohd, dash-dotted and dashed hues are for the results with different treatments of the 
mass-sheet degeneracy, specifically, with true cluster mass distribution, with < k >— in 
the ring at 5Rc, and with < k >— within 5Rc, respectively. The cluster parameters are 
the same as those in Figure 2. 
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Fig. 7. — The height distribution of the main-cluster-peak measured from noisy weak-lensing 
convergence maps. The upper and the lower panels are for the isothermal and the NFW 
clusters, respectively. The solid line in each panel is for the result calculated from Eq. (31). 
The histogram is for the result from our Monte Carlo simulations. The vertical dotted line 
indicates the value of = 5.6, the true peak height of the considered cluster. The cluster 
parameters are the same as those in Figure 2. 
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Fig. 8. — The spatial offset distribution of the measured main-cluster-peak position to the 
true center of the cluster. The sohd and dash-dotted hues are for the results calculated from 
Eq. (34) for the isothermal cluster and the NFW cluster, respectively. The diamond and 
plus symbols are for the corresponding results from Monte Carlo simulations. The cluster 
parameters are the same as those in Figure 2. 



-33- 




Fig. 9. — Comparisons for three NFW clusters of about the same mass but with different 
density profiles. The lines are the results from our theoretical calculations. The histograms 
and the symbols are the corresponding results from Monte Carlo simulations. The solid line 
in each panel is for the cluster with Rg = and central K/aQ = 3 x 5.6. The dashed line is 
for the cluster with Rg — 1.8i?c and central K/gq — 1.8x5.6. The dash-dotted line is for the 
cluster with Rg — 5Rc and central K/ao — 5.6. Fig. 9a (upper left panel) shows the profiles 
of K/ao for the three clusters. Fig. 9b (upper right panel) is for the height distributions of 
the main-cluster-peak for the three clusters, where the vertical dotted lines mark the values 
of the true peak height of the three clusters, and the vertical dash-dotted lines show the u 
value where p(z/) is peaked for each of the three clusters. Fig. 9c shows the spatial offset 
distributions of the main-cluster-peak position for the three clusters. 




Fig. 10. — The differential peak tieiglit distribution in 1 deg^. Tlie liistograms are for tlie 
results averaged over 16 ray-tracing simulations with 3x3 deg^ each by White and Vale 
(2004). The lower and upper histograms correspond to the results without and with noise 
included. The dotted line is for the theoretical result considering only NFW halos without 
noise [Eq. (40)]. The long-dashed line is for the result of dark matter halos including the 
noise effects on the their peak heights calculated from Eq. (41) and Eq. (31). The solid line 
is for the prediction from Eq. (35)-Eq. (38) including the full noise effects. 



